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Abstract 

The hadronic decay of radially and orbitally excited charmonium above charm 
threshold by 3 Po pair creation and chromoelectric flux-tube breaking is discussed in 
an harmonic oscillator approximation. We find independent evidence from a study 
of widths for a 2S admixture in the predominantly ID state ^(3770), and explore 
the possibility of metastable radially excited 2 3 -Fn,i,2 states being a source of the 
anomalously large production of tp at the Tevatron. At least one of them is expected 
to be narrow as a consequence of the existence of nodes in the radial wave function. 
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1 Introduction 



Recently ip enhancement at CDF [|IJ has generated considerable interest. This may cur- 
rently be the largest discrepancy between the predictions of the standard model and ex- 
periment, and arises because the behaviour of QCD in the strongly interacting region is 
inadequately understood theoretically. There are four possible sources of enhanced ip dis- 
cussed in the literature 0, |, g, || : 2~+ cc below the D*D threshold @, |] (D = 0~ and 
D* = radial ^-states (i.e. 2P charmonia) []3], Q, hybrid charmonium || and ip with a 
ccgg-component ||. Investigation indicated that 2~ + is an unlikely candidate @, |3|] unless 
ip contains significant 3 Di in its wave function. Detailed study of fusion and fragmen- 
tation contributions to the production of radial x~states indicated that they could explain 
the ip excess if they are very narrow || || (~ 1 MeV). We investigate whether this is likely 
or even possible in a QCD inspired model known to be successful for light hadrons ||, 
namely the non-relativistic flux-tube model of Isgur and Paton [[/J , which is related closely 
to the phenomenologically successful [0] 3 Po~ m odel ||. The extension of earlier work on 
light mesons |j] to charmed hadrons is novel in its own right, and will form the basis for a 
subsequent study [TI| of hybrid charmonium. 

The strategy is to study first the DD width of -0(3770). This is the cleanest example ITT 



being a natural 1D-2S candidate. The only real uncertainty is the lD/2S-mixing, though 
e + e~- and 7-transitions suggest a dominantly ID state. We find that the width fits with 
parameters determined elsewhere || [TJ|, O, 111, |15j; and discover that the next potential cc 



state of higher mass, the 3S [ 16fl , would not enable sensible fitting of the -0(3770) widths. 

It is encouraging that the width of -0(3770) fits well in this picture. In fact, it fits in so 
well that we find the need to constrain the model parameters by also fitting energetically 
higher lying states, such as -0(4040), 0(4160) and ■0(4415), to experiment. 

When the model is applied to calculations of the widths of radial x~states the results 
depend rather critically on the mass of these states, and hence the phase space available for 
their DD and D*D decays. The 2 ++ decays in D-wave into DD, D*D, while 1++ -> DD* 
are likely to be near to threshold. As anticipated in ref. || we find the 2 ++ generally quite 
narrow (~ 0.5 — 5 MeV) in the range of predicted masses. The 1 ++ and ++ widths are 
very sensitive to the nodes of the radial \ wave function. It is possible that they are very 
narrow (~ 1 MeV), but this would be a coincidental conspiracy, which nevertheless can be 
attained near the nodes. We find that for sensible parameter solutions to ^(3770), it is 
likely that some of the (0, 1,2) ++ widths are considerably reduced, which is sufficient for 
the enhanced production of ip at CDF. 

The basic structure of the paper is as follows. 

In §2 we discuss the charmonium system and introduce the flux-tube model and the 
parameter values used. In §3 we fit the ■0(3770) to experiment and compare the range of 
possible mixing angles with those obtained from e + e~ annihilation. We are encouraged by 
the consistency of these results, which prompts us to look at the higher radially excited 
states in §4, and observe to which degree their experimental signatures can be accounted for. 
The parameters obtained are used in §5 to look at radial x~states and to examine whether 
they may be a source of the anomalously large production of ip at CDF. We conclude by 
outlining the experimental consequences. 
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2 Outline of cc phenomenology 



Earlier work [11], [lo[ [17], [L8| on the charmonium system JT9], [Hj, |2T[ is extended by taking into 
account flux-tube breaking || |22|. We perform an analytical calculation of radially and 
orbitally excited charmonium decay amplitudes to DD, D*D, D*D*, D S D S , D*D S , D*D* 
and D**D (D** = (0, 1, 2) + ) in a S.H.O. wave function approximation. This approximation 
allows comparison with previous work || [15|, |17], [18| in the 3 Po-model, improves the ability 



to handle future parameter changes, and is known to be an excellent for charmonia and 
charmed mesonsQ The calculation is performed in the approximation where the inverse 
radii j3 of the outgoing mesons are identical. The main reason for this simplification is 
that the D** , D s , D*, D* and D are expected |J to have similar /3's, and that the flux-tube 
and 3 Po-models correspond closely in this case (see the next paragraph). Results do not 
depend critically on this simplification, making it unnecessary to concern ourselves with 
small corrections. 

The flux-tube model decay amplitude || by pair creation differs from the familiar 3 Po- 
model amplitude, although they coincide in the case when (5 is a universal constant for all 
outgoing mesons. This is shown in Appendix A, § [A.2| . In this work, we allow variation of 
(3. Hence the two models are not identical, though their qualitative features are similar [|J, 

!• 



17, m 



allowing comparison with earlier studies [y, [T 

The only free parameter in the model is the overall normalization of decays g^|^oo 
(eqn. [IT], Appendix B). This dimensionless factor is common with light meson decays and 
was phenomenologically found to equal 0.64 M for creating light quark pairs with u, d and 
s flavours^} We adopt it here; thus implying that cc widths are predicted independently. 

From the ISGW non-relativistic fit [[^J to spin-averaged meson masses we take the string 
tension b = 0.18 GeV 2 , and the constituent-quark masses m Ui d = 0.33 GeV, m s = 0.55 GeV 
and m c = 1.82 GeV. Meson masses are taken to be the PDG masses [24 , and where not 
available (which is also the case for 3 P 1 / 1 P 1 mixing angles) motivated by spectroscopy 
predictions [pT6fl , adjusted in absolute value relative to known masses. The procedure for 

and |25], §2]Q. 



the calculation of widths is outlined in refs. 



The experimental total decay widths in MeV are [24 



VK3770) 
0>(4O4O) 
V>(4160) 
0>(4415) 



DD t 

DD t,D*D t,D*D* t,D s D s 
DD, D*D, D*D*, D S D S , D*D S 
DD,D*D,D*D*,D a D« 



D* S D S 



D*D*,D**D 



23.6 ± 2.7 

52.0 ± 10.0 

78.0 ± 20.0 

43 ± 15 



(1) 



where we indicated the dominant^ decay modes considered in this work, and f indicates 
modes that have been observed ||24|| . The modes indicated in eqn. [I] are assumed to 
dominate over decays to cc + mesons. In the next section we shall calculate these 



1 When expanding full numerical wave functions the overlap with the corresponding radially excited 
S.H.O. wave function is typically found to be dominant, e.g. 99.8% (D, D*), 97.7% (V>(404Q)) and 97.0% 
(-0(4415)) Jl5| for appropriately chosen j3. The nodal positions in wave functions and decay amplitudes are 
also well approximated fllq l . 

2 Suppressions of 1 JlT] , Q, m u /m s |l7| ] and (m u /m s ) 2 [ pT"| in amplitude have been employed for ss 
creation relative to uu, dd creation. We adopt the first. 

3 We employ for the masses M [| the values 1.97 GeV (D,D*), 2.07 GeV (D S ,D*), 2.44 GeV (£>**), 
3.76 GeV (^(3770)), 4.015 GeV (0(4040)), 4.15 GeV (0(4160)), 4.39 GeV (0(4415)) 

4 OZI forbidden decays are are expected to be small (e.g. experimental indicattions are that 0(3770) — * 
07T7T has a width of - 20 - 80 keV |7)). Also, 0(4415) decays to e.g. D(2S)D, D**D*, D*(2S)D* and 
D**D* are expected H |3 to be kinematically forbidden. 
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Table 1: The outgoing meson on-shell CM momentum ps in GeV. 



^(3770) -> DD 



0(4040) 
0(4040) 
0(4040) 
0>(4O4O) 



0(4160) 
^(4160) 
^(4160) 
0(4160) 
0(4160) 



DD 

D*D 

D*D* 

D S D X 



DD 

D*D 

D*D* 

D S D S 

d:d„ 



Pb 



0.26 



0.77 
0.57 
0.21 
0.45 



0.92 
0.75 
0.54 
0.67 
0.41 



V>(4415) 
^(4415) 
^(4415) 
^(4415) 
^(4415) 
^(4415) 
^(4415) 
^(4415) 
0(4415) 



DD 

D*D 

D*D* 

D S D S 

D* S D S 

D* S D* S 

D*X+D 

D\% L D 

D{% H D 



Pb 



1.18 
1.06 
0.92 
1.00 
0.84 
0.65 
0.44 
0.52 
0.46 



contributions to the hadronic width in our model, and sum them up in order to fit them to 
the total experimental widths in eqn. |l[ 



3 Fitting the overall widths to experiment 

3.1 ^(3770) 

For some time there has been a need to obtain an overall picture of excited charmonium with 
modern parameters. The process 0(3770) — ► DD produces a D-meson of momentum pb = 
0.26 GeV (see table |l|). The decay happens far from where the widths vanish, i.e. from the 
nodes in the lD/2S-amplitude, for all phenomenologically relevant (/3a, [3). This, the fact 
there is only one hadronic decay mode, and that the width has been estimated successfully 



before |TT| , 26], make it the cleanest charmonium candidate above charm threshold. We 
shall now proceed to discuss this state in more detail. 

The J/ip has unambigously been identified as dominantly IS, so the cc state ■0(3770) 
must be a higher radial or orbital excitation: most probably, on mass alone, a 1D-2S 



mixture |IE |. This can independently and non-trivially be established by considering its 
width. We find that the -0(3770) width is consistent with it being ID or 2S, i.e. with 
a 1D-2S mixture. The alternative "gedanken" assumption that 0(3770) is a 3S or 4S 
state leads to a contradiction: this would imply a width to DD too small to be consistent 
with experiment (for PaiP ~ 0.1 — 1 GeV). Hence we believe that the 1D-2S nature of 
0(3770) established from width considerations, being consistent with mass spectroscopy, is 
significant. 

Noting that the ip (3685) is predominantly 2S, we expect 0(3770) to be mainly ID, with 
a 2S admixture^ as |0(377O) ) = cos9 \1D) + sm9 \2S). In figure 1 we indicate regions in 
parameter space where 0(3770) can be fitted to experiment as a 1D-2S mixture for various 
9. The fact that the width can be fitted at all is non-trivial, as will be briefly outlined. 
If we mix the ID and 2S decay amplitudes as a new amplitude / = cos 9 ID + sin 9 2S, 



5 We neglect IS admixture for the following reason : Although our calculation shows that taking ^(3770) 
as pure IS is sensible for 0a ~ 0.1 — 0.5 GeV and ~ 0.2 — 0.8 GeV, and it is known that IS mixing 
in 0(3770) can be significant |]36f for spin-dependent forces; it is small in coupled channel treatments 
p[ p9| , p7| . The latter models successfully reproduce the experimental 2S admixture. 



5 



then f 2 < (ID) 2 + (2S) 2 ; with the maximum of |/| occuring at tan# = 2S/1D and the 
minimum / = at tan 9 = —1D/2S. Clearly there is a large range of widths that the 
mixed state can have as 9 is varied. There is still a restriction, though: in regions where 
(ID) 2 + (2S) 2 < 24 ± 3 MeV, the 0(3770) total width f 2 < (ID) 2 + (2S) 2 cannot be fitted 
for any 9. However, there is a large region where / is large enough to reproduce experiment. 

We can restrict the parameter space (Pa, P) by incorporating theoretical prejudice about 
the /3's. D, D* and D** have p in the region^ 0.30 - 0.70 GeV; and ip(2S) and 0>(1L>) have 
Pa = 0.30 — 0.55 GeV. From figure 1 we note that '0(3770) can be fitted for any 9 in this 
region of parameter space, although a sizable 2S admixture with sin 9 ~ 30 — 45 % appears 
to be preferred. This will be reinforced when we fit more massive states in 



33.2 



It can be estimated^ from the observed e + e -widths of ip' and '0(3770) that either sin9 
= 17% [g7] or 23% or sin 9 = -49% [0] or sin# = -44% J28|, employing our convention 
for 9. Clearly our analysis prefers sin# > 0, concurring with findings on the detailed 
analysis of ip-Kir |27|| and cc radiative [28| decays. This also in agreement with theoretical 
calculations in the coupled channel formalism, which find sin9 = 17% or 20% [29 . 

Given that the overall normalization of decays was taken from light meson spectroscopy 
as part of a unified study, in contrast to earlier cc calculations |L5], [L7|, [L8|, ^2], |30[, we do 
not regard the disagreement of our estimate with that from e + e~-widths as significant. 
Our determination of the 2S admixture in 0(3770) is hence broadly consistent with its 
e + e~ width and the phenomenology of the ip . We thus consider (Pa, P) consistent with 
a realistic and dominant ID component in 0(3770), but with a sizable 2S admixture of 
~ 30% in amplitude. 

In summary we conclude that this model of decay by flux-tube breaking can successfully 
account for the 0(3770) as a 1D-2S mixture in this study of widths, consistent with both the 
a priori theoretical estimates for (Pa, P) and the restrictions on 9 from e + e~-widths. These 
conclusions based on its width are also in line with mass spectroscopy. In the next section 
we shall find that higher states can also be described with varying degrees of success. 



3.2 Higher mass charmonium 

As far as states more massive than 0(3770) are concerned, we are led to identify them as 3S- 
2D mixtures, or higher radial excitations, due to the existence of convincing candidates for 
IS, 2S and ID states. The 0(4160) and -0(4040) (which may be either simple resonances 
or a collection of states), should at least contain the 3S and 2D states, otherwise the 
0(4415) must be 3S or 2D, leading to a contradiction: the IS, 2S, 3S, ID and 2D wave 
functions produce 0(4415) total widths too large to be consistent with experiment (for 
Pa, P ~ 0.2 — 0.8). Hence by default the natural interpretation of 0(4160) and -0(4040) 
is that they are 3S-2D mixtures. They can be fitted consistently with experiment in the 
phenomenologically relevant region Pa, P ~ 0.3 — 0.6. Also, 0(4415) can be fitted as 4S 
for Pa ~ 0.3 — 0.6 GeV and P ~ 0.4 — 0.6 GeV. We shall hence adoptf] these assignments, 



which are in agreement with mass spectroscopy [16, 31 



6 It is estimated that (3 D . D * = 0.39 GeV ||, 0.40 GeV JTJ, 0.42 GeV @, 0.49 - 0.55 GeV @ or 0.54 
- 0.66 GeV §; and p D „ = 0.34 GeV @, 0.42 - 0.49 GeV or 0.45 - 0.54 GeV § (on average 0.07 GeV 
lower than p DtD .). For cc the estimates are : p is = 0.66 GeV g% 0.75 GeV [[TJ or 0.57 GeV @; Pi P = 
0.50 GeV |§;'/3id = 0.37 GeV jfj, 0.45 GeV @ or 0.48 GeV|2| 



The experimental uncertainties in e + e~-widths |24 translates into ~ 4% uncertainty in sin#. 
8 We assume no mixing for higher charmonium states for simplicity : it is generally considered to be 
small pH, e.g. ~ 10% pil in coupled channel models. 
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Table 2: A selection of widths in MeV fitting experiment (see §[3]). We assume a common f3 
(in GeV) in order to reduce parameters. We allow 0(3770) to have P^- 3770 ' = Pa + 0.1 GeV 
(see § [3.2|) , and a 1D-2S mixing angle 9, in the decay -0(3770) — > DD. We indicate 0(4160) 
as either 3S or 2D, and the widths of the various decay modes 0(4160) — > DD, D*D, D*D* 
(P- and F-wave), D S D S ,D*D S . 







0(3770) 




0(4160) 


D*D* 








Pa 


P 


sin 9 


Total 




DD 


D*D 


P 


F 


D S D S 


D*D S 


Total 


0.50 


0.54 


0.3 


23 


2D 


32 


1 


6 


20 


4 


11 


74 


0.50 


0.54 


0.3 


23 


3S 


13 


4 


34 





1 


20 


72 


0.42 


0.52 


0.4 


24 


3S 


17 


44 


3 





1 


2 


67 


0.30 


0.39 


0.5 


21 


3S 


1 


22 


52 





5 


1 


81 



Table 3: As in table | but for 0(4415) decaying to DD, D*D, D*D* (P-wave), D S D S , 
D*D S , D*D* (P-wave), D* 2 * ++ , D\% L (S- and D-wave), D*%h (S- and D-wave). Here L,H 
indicates the low and high J p = 1 + mass states respectively, with a 3 -Pi/ 1 Pi -mixing angle 



of 9 = 


-41° 


lfj (see eqn. 


13], Appendix B). We take the P 


IXh mass to be 2.45 GeV. 


Pa 


P 


DD 


D*D 


D*D* 


P S P S 


D*D S 






D 


** 
l+£ 


^l+H 


Total 




















S 


D 


s 


D 




0.50 


0.54 


1 


3 


1 


1 


2 





14 





17 


1 


1 


41 


0.42 


0.52 


8 


22 


11 


1 








1 





1 


1 





45 



We expect the -0(4160) to be the easiest to describe. This is because there is an experi- 
mental constraint on the width ratios of 0(4040) |T7| (see below), casting doubt on whether 
it is conventional cc. Also, the '0(4415) has numerous decays modes (including decays to 
D**D), making its narrowness |L7| potentially difficult to understand. 

It is theoretically expected that the magnitudes of Pa for higher excited states are 
smaller than those of lower states (e.g. 0^(3770) in this case) by ~ 0.1 GeV ||. Allowing for 
this, we findQ that all three of the higher states can be accommodated with a 2S admixture 
in 0>(377O) of ~ 30 — 45 % in amplitude. It can be seen by comparing figures 1 and 2. 

If we choose parameters consistent with a sizable 2S component in ^(3770), figure 2 
indicates that it is possible to fit ■0(4160) as 3S, and to a lesser extent as 2D, in a wide 



range of theoretically acceptable (Pa,P)- If we require P near to the ISGW value p3 
(Pa, P) = (0.30,0.39) GeV gives a good fit for most states. Some of these values also allow 
fits of ^(4415) (see tables | and |), e.g. (p A , P) = (0.42,0.52) GeV, where we find excellent 
agreement. The point (Pa, P) = (0.50,0.54) GeV, which fits -0(4040) in addition (see below), 
also lies in this region. 

We now proceed to discuss the -0(4040). The issue of whether it can be understood 
as ordinary charmonium has historically been controversial, due to its anomalously large 
D*D* branching ratio. The large experimental ratios p4] Ri = T(D*D) / T(DD) = 5 - 



20 and R 2 = T(D*D*) / T(D*D) = 1-2.5 can conceivably be explained by looking for 
nodes (zero's) [|TTJ] in the decay amplitudes corresponding to decays into DD and D*D. 
This can arise if the kinematics of the decay are controlled by nodes arising from the 



9 This result holds for 0(4040) as 3S, and to a lesser extent as 2D; for ^(4160) as 3S, and to a lesser 
extent as 2D; and for a 4S 0(4415). 
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Table 4: Total widths in MeV at specific {(3a, 0) (in GeV) providing an excellent fit to 
^(4040) (see M). -0(3770) is fitted with a 1D-2S mixing angle 9 assuming /^ (3770) = 



(3a + 0.1 GeV. The -0(4040) width ratios i?i ;2 defined in § )3~2] are indicated along with the 
widths of the various decay modes. ■0(4040) is constained to be 3S, implying that -0(4160) 
should be 2D. All decays are in P-wave. 











0(4040) 








0(3770) 


V>(4160) 


(3a 




DD 


D*D 


D*D* 


D S D S 


Total 


Ri 


R2 


sin# 


Total 


Total 


0.50 


0.54 


2 


15 


23 


5 


45 


9 


1.5 


0.3 


23 


74 


0.24 


0.23 


4 


26 


34 


2 


66 


6 


1.3 


0.1 


22 


41 


0.50 


0.20 


1 


14 


29 


8 


52 


15 


2.0 









radial wave functions of the 3S or 2D states. There is thus the possibility of the DD and 
D*D amplitudes being near enough to two[^] different nodes or to the same node [IT]. We 



proceed to seek cases {(3a, (3) where Ri and _R 2 are consistent with experiment. Although 
R lj2 can be fitted for 0(4040) as 2D with {(3 A ,(3) ~ (0.3,0.2), the total width is too small, 
leading us to the conclusion that a consistent picture of -0(4040) can only be obtained with 
it as 3S. Our search provides three areas in parameter space where the -0(4040) may be 
realized as 3S (see table f| and figure 2). If in addition we require the -0(3770) and 0(4160) 
to be consistent with experiment, -0(4160) to be 2D (as a consequence of -0(4040) being 
3S), and the parameters to be consistent with theoretical estimates, we are restricted to 
{(3a, (3) around (0.50,0.54) GeV. This region represents the only area where we can fit all 
experimental data on excited charmonium decays, including the -0(4160) and 0(4415), with 
the 0(3770) a 1D-2S mixture as usual (see tables ^| and (|). The region also appears to be 
related to the one found by Le Yaouanc et al. ]T7|], who also identified^ -0(4040) as 3S. This 
implies that it is possible to fit -0(4040) consistent with experiment without contradicting 
the phenomenology of the other 0-states, although it severely constrains the parameters, 



requiring considerable coincidence of parameters [0 and fine tuning of to ~ 1 — 3%, 
in order to obtain the correct branching ratios. 

The fits obtained for {(3a, (3) = (0.42,0.52) and (0.50,0.54) GeV above are remarkable in 
the sense that we can accurately fit three amplitudes (i.e. 0(3770), 0(4160) and -0(4415)) to 
experiment with two parameters (3a and (3, with the additional conclusion that -0(4040) can 
possibly be fitted, and furthermore that the assigned nS, nD classifications are consistent 
with mass spectroscopy. 

In summary we conclude that the higher states may also be accommodated in this 
model. We have thus established another rung on the cc ladder, highlighting the difficulties 
associated with -0(4040) and ^(4415). 



10 We showed that there exits no (Pa,P) for which for which the momenta ps — 0.77 GeV (DD) and 
Pb — 0.57 GeV (D*D) lie near different nodes in the decay amplitude. 

11 They obtained (3a — P — 0.44 GeV, but with old parameter values and a mass formula neglecting the 
difference between m c and m u . 
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4 2P charmonia and the ip anomaly at CDF 



The experimental production rate of ip' is 30 times larger than theoretical estimates [p]] at 
large transverse momentum in pp collisions at the Tevatron. One proposed source H || is 
the production of 2P charmonia followed by their radiative decay 2 2S+1 Pj — » 2 3 S\ + 7, 
enhancing the ip rate. This can only happen if the 2P charmonia have small hadronic 
widths. These states are expected [[L6|, |31| to lie above the PP-threshold (3.73 GeV), 
probably above the P*P-threshold (3.87 GeV), and possibly above the P s P s -threshold 
(3.94 GeV) but beneath the P*P*-threshold (4.02 GeV). In these circumstances we find 
that their hadronic widths may be small. 

If the 2P states lie below the P*P-threshold (say at 3.85 GeV), then the 2 3 P 2 and 
2 3,1 Pi are narrow : 2 3 P 2 (because of D-wave phase space) has width 0-4 MeV, and 2 3,1 Pi 
is narrow because only decays to D*D are allowed. 

If the 2P states lie above the P*P-threshold the 2 3 P 2 has the unique property that the 
total width remains small (because of D-wave phase space) at - 7 MeV (3.90 GeV) and 1 
- 14 MeV (3.95 GeV). Figure 3 indicates that in the theoretically reasonable region where 
(Pa,P) ~ (0.3,0.4) - (0.4,0.6) GeV these widths can be £ 1 MeV if (p A ,P) is tuned to 
only ~ (15,20) %. This is of special interest for the resolution of the CDF anomaly. The 
2 3,1 Pi and 2 3 Po states only have small widths if the kinematics of the decay cause the 
amplitude to be controlled by the nodes arising from the radial wave functions of the 2P 
states (otherwise the amplitude can be substantial, even near thresholds). 

• For decays of 2 3 ' 1 P 1 the S-wave nodes occur in a line of parameter space (Pa,P) ~ 
(0.3, 0.4) — (0.45, 0.6) GeV coinciding with narrow 2 3 P 2 (see figure 3). It is fascinating 
that this region coincides with D-wave nodes, reinforcing the narrowness of 2 3,1 Pi. 

• For decays of 2 3 P the S-wave nodal line occurs for (Pa, P) ~ (0.4,0.3) — (0.5,0.6) 



GeV (see figure 3). Estimates |T6[ suggest this to be the lightest 2P state, and hence 
possibly narrow (i.e. below the P s P s -threshold, which ensures that 2 3 P — > D S D S 
does not overshadow the decay). 

Hence we expect at least one of 2 3 P 2 , 2 3 Pi or 2 3 P to be narrow; the 2 3 P 2 being 
most likely, since obtaining 2 3 P and 2 3 Pi widths of < 1 MeV requires fine tuning (Pa, P) 
respectively to ~ (4,7) % and ~ (1,1) % in figure 3. Theoretically we expext /3 2 p ~ 0.3 — 0.45 
GeV, tantilizingly in accord with the values where we expect nodes to occur. 

Both the best regions in fitting experiment (i.e. near (Pa, P) = (0.42,0.52) and 
(0.50, 0.54) GeV) can be consistent with one of the 2P states being very narrow (~ 1 
MeV), depending sensitively on their masses. 

If an amplitude lies near a node, it becomes especially sensitive to differential amounts 
of phase space e.g. for D°D° and D + D~ final states. This sensitivity may be decreased 
if other small decay modes far from a node swamp the decay. Since dps ~ 1/pb, the 
change in available phase space is largest when the decay is nearest to threshold, so that 
amplitude variations^ are pronounced near thresholds. Hence, when 2P masses are known 
experimentally, these effects should be incorporated near the various thresholds if there is 
nodal suppression of an amplitude. 



12 We have not considered these effects in §3.1. They are expected to be largest for ^(3770) — > DD, 
because it is near threshold (see table [l]), and the dominant hadronic decay mode. However for realistic 
(Pa, 13) it is far from a node, with typically 5 - 15 % variation in width due to varying the phase space of 
decays to D°D° and D + D~. So the effect is small compared to variations in ((3a, (3). 
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5 Conclusions 



Our results from hadronic decays of charmonium are as follows : 

• 0(3770) — > DD can be understood with the ■0(3770) as a mixture of 2 3 S' 1 and 1 3 Di. 



We find that the width fits with parameters determined elsewhere || [12|, y| [14 
and discover that the next state of higher mass, the 3S |LJJ, would not enable sensible 
fitting. We favour a sizable 2S admixture in 0(3770) (~ 30% in amplitude), consistent 
with e + e~-annihilation. This is also in observational agreement with higher states. 

We find further encouragement that the widths of 0(4160) and 0(4040) fit well if they 
are 3S-2D mixtures. We note that they do not admit solutions as 2S-1D states. 

0(4160) -> DD,D*D,D*D*,D S D S ,D*D S can be understood with the 0(4160) as 
either 3 3 S\ or 2 3 Z?i, depending on the preferred value of (/3a, (3) in §§. Study of the 
branching ratios of this state at a r-charm factory could be a rather critical test of 
its wave function composition. 

0(4040) cannot be understood as 2 3 D 1 , due to an inability to reproduce simultane- 
ously the total width and correct branching ratios into DD,D*D,D*D*. These are 
experimentally (neglecting statistical uncertainties) 1 : 8 : 14. Remarkably, the state 
can be understood as 3 3 Si in three regions of parameter space (see table |j). This, 
however, severely constrains parameters. Previous interpretations of 0(4040) as cc 



15] may hence now be invalid due to a change in preferred parameters. This leaves 



open the possibility of 0(4040) being a D*D* molecule |19| or a threshold effect p2 | 



(the state is 10 - 30 MeV above the D*D* threshold with a width of 40 - 60 MeV, 
and can hence not be treated well in the narrow resonance approximation). There 
may be several 1 states in the 4.0 - 4.3 GeV region, due to the possible additional 



existence of a hybrid meson of mass ~ 4.2 GeV |3T], |33f, but whose e + e production 



is suppressed due to the radial wave function vanishing at the origin [ |3T|j . 

The 0(4415) is pinned down to having a 4S wave function. Its total width tends 
to be small enough only in restricted regions, suggesting that the narrowness [|T^] of 
0(4415) remains an interesting issuef^]. It is not easy |30]] to identify dominant decay 



modes without fixing ((3a, (3) and the degree of mixing, as was effectively done in the 
literature [17|, [18], [1^], in contrast to earlier expectations |17|, [18| where D^X+D was 
indicated as possibly significant. 

Where energetically allowed, decays to D s , D* are suppressed by a flavour factor of two, 
relative to decays to D, D*. This, together with P-wave phase space, conspires to make the 
branching ratios of 0(4040), 0(4160) and 0(4415) to D and D* consistently larger than 
those to D s and D*. This is violated for 0(4040) — > D S D S , where the lack of a corresponding 
experimental width estimate may indicate additional suppression of ss creation [11, T7J (see 
table £| and footnote ||). 

Radial x~states may be the source of a significant cross-section at the Tevatron if 
their widths are narrow. Their widths depend rather critically on their masses, and hence 



13 The narrowness can be increased by increasing the radial excitation. Studies of the 0(4415) width 
[fl7| , [l8| , [2^] and mass (^9) tend to find a 4S assignment. A 5S assignment (with an accompanying 4S 
■0(4160)) has, however, been suggested This possibility is found to be consistent with experimental 

widths by our study, and can as such not be ruled out. 
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the phase space available for their DD and D*D decays. The 2 ++ decays in D-wave into 
DD, D*D, while 1 ++ decays to DD* are likely to be near to threshold. As anticipated in 
ref. we find the 2 ++ generally quite narrow (~ 0.5 — 5 MeV) in the range of predicted 
masses, with the narrower widths corresponding to D-wave nodal suppression. The 1 ++ 
and ++ widths are very sensitive to the nodes of the radial x wave function. It is possible 
that they are very narrow (~ 1 MeV), but this would be a coincidental conspiracy. For 
++ below the /^.D,; -threshold, nodal suppression is quite consistent with ^(3770) decays. 
We thus find that for sensible parameter solutions to ?/>(3770), it is likely that some of 
the (0, 1,2) ++ widths are considerably reduced. This would suggest placing emphasis on 
searching for the 2 3 P 2 in events containing ip . We predict that if its hadronic width is 
found to be ~ 1 MeV, then the 2 3 > 1 P 1 may also be found to narrow (depending sensitively 
on their mass). 

The overall consistency of excited charmonium with observation vindicates the expec- 
tation that the pair creation amplitude is similar for charmonium and light meson decays 
(in each case light quarks are created). In view of the slightly exaggerated 2S admixture we 
obtain in -^(3770), this amplitude may have to be adjusted slightly. This procedure would 
however not affect the nodes in the \ decay amplitude, leaving the (/3a, P) where x~states 
are narrow unaltered. 

We highlight the need for improved data on excited charmonia : to confirm whether they 
are true resonances, to study their branching ratios into various channels and to discover 
whether some may be narrow. Rather clear signals, in particular in e + e~-annihilation, may 
also reveal hybrids and molecules. 
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A Appendix : The coincidence of flux-tube and 3 Po — 
model decay amplitudes in limiting cases 

For meson decays to mesons the 3 Pq and the flux-tube models coincide in the limit of an 
infinitely thick flux-tube, i.e. where the pair creation amplitude is constant all over space. 
This is indicated explicitly in §[A. 1|. It is possible to make a stronger statement : The 3 Pq 



and flux-tube models coincide when (3 is a universal constant for all outgoing mesons, even 
if the flux-tube has finite thickness. This is demonstrated in § [A.2| . 

We focus on two initial quarks of mass M, with pair creation of quarks of mass m. 



The decay amplitude of an initial meson A into final mesons B and C can be shown p5 
Appendix A] in the rest frame of A (where = 0) to be given by 



ac ■ - ^3^3 



^(2^) V(p B + Pc ) -TriA 1 BC) fi avour Tr(A Bcr 1 C) spm - 

x f d 3 r A d 3 y ip A {*A) exp(i — ^-ttPb ■ r A )^(r A , y ± ) 
J m + M 

x (*V rB + *V rc + -—^p B )^ B (r B )r c (rc) + (B - C) (2) 
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Here (B <->• C) indicates a term obtained by interchanging the flavour and spin matrices 
B <-» C and momenta p B *-> pc in the first term in eqn. |2|. For cc decaying by the creation 
of a mm, (id or ss pair, only one of the terms contribute (due to Tr(A T BC) fi avour ); and the 
space part of the term differs by a sign from the space part of the displayed term in eqn. 
0. So from now on it is sufficient just to consider the the displayed term. 

The pair creation position y is measured relative to the CM of the initial quarks, and 
y_i_ = — (y x ta) x is a perpendicular "component" of y to the initial QQ-axis ta- The 
Qq - axes of the final states B and C are i\b = + y and rc = i"a/2 — y respectively pi. 

A.l Spatially constant pair creation 

The pair creation amplitude (or flux-tube overlap) is 




The thickness of the flux-tube is related inversely to /. 
quantities and the structure of eqn. |3] may be found in ref. 
The estimated values / = 1.1 and Aq = 1.0 [21 
thick flux-tube with / = corresponds to the 



(3) 



A detailed discussion of these 
A21] and ref. |3 



eqn. 



34J are used in this work. The infinitely 
In this case 7(1% y±) 



Po~ m odel |]Tyj) . in tms case 7(1% yjj is 
a constant (which should be normalized to be non-zero), and the gg-pair is created with 
uniform amplitude anywhere in space. Fourier transforming (a relevant part of) eqn. |2| 
yields 



— (27r) 3 5 3 (p B + p c ) e a - / d 3 r A d 3 yipA(r A )exp( y i — — p B • r A ) 

2 J m + M 

2m 

X (lV rB + lV rc + ^J^Pb) ^bWc^c) = 

-i5'\p B + p c ) e a - I rf 3 k A (k A + p B )^ A (k A )^(k A + — - p B )i)* c (k A + — — p B ) (4) 
J m + M m + M 

where we defined ipA,B,c{ r ) = (2vr) -3 / d 3 k exp(ik ■ r)ipA,B,cQ<-)- 

Taking the 3 P ~ m °del amplitude |35], p. 120], with the gg-pair created with equal but 
opposite momenta 

TVT / ^i^A^i^i + fe-P-i) i 3 (k, + ka-p B ) * 3 (k 2 + k4-p ) i 3 (k 3 + k 4 ) 

x n„(k 3 - k 4) M±i±) m^ip-) (5) 

2 m + M m + M 

we can see that it equals the last line of eqn. f| jy| (with pa = 0) when we make a change 
of variables to k A = (ki — k 2 )/2 and kcM = ki + k 2 . Here kYi a (k) = e CT ■ k was used. 
We have thus connected the flux-tube and 3 Po~ m °del notation explicitly; showing that the 
two models coincide for an infinitely thick flux-tube; and verifying the expectation that a 
zero momentum qq-pair corresponds to complete freedom in creating the pair anywhere in 
space. 
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A. 2 Outgoing mesons of equal "size" 

The incoming and outgoing mesons have S.H.O. wave functions with inverse radii f3 A and 
(3 respectively. We assume the outgoing {D, D* , D s , D*} and D** wave functions to have 
angular momentum L=0 and L=l respectively. In the lowest radial state the wave functions 
are: 

Mr) = Mlt l Y LMl {v) exp(-/?£r 2 /2) M L = -^{1, ^ -(3 L , -^=0*} (6) 

All wave functions are properly normalized (with the brackets referring to S, P and D states 
respectively). Performing the differentiation in eqn. |2] gives for S + S final states 

e* • (*V rs + iV rc + m + M PB) ip* s ( r B)ip*s( r c) = 

Ml eM-P 2 C-f + y 2 )) HP*t a + -^j-fvs) (7) 

4 m + M 

and now proceeding to perform the y-integration gives 



// to Ztti 
d 3 yi(r A , y ± ) x (eqn. 7) = A° m \ h^Afl f {-i(3 2 v A + — -p fl ) exp( 
V 7r m + M 



2„2 

A ' 



m + M 

-il _ If, , «2 2 2 

7 = d y exp(-/3 y - — y ± 



+ /6/2) 



Similarly for P + S final states 



2772 

(*V rs + 2V rc + — -p B ) Vp(r B )^s(rc) = A/^Vp exp(-/3 2 (^ + y 2 )) 

m + M 4 

x [i V^ e k fl + (yW^) + y^(y)) HAa + J^mPb)) (9) 

and performing the y-integration, terms linear in y (in eqn. |9|) vanish 

J d 3 y-f(v A ,y ± ) x (egn.fj) = 

A^UpUs 7° [*'/J^M f + y W^) + J^tjPb)} ex P (-^) (10) 

Since the 3 P ~ m °del corresponds to taking / = in the amplitude above, we see 
that the flux-tube model amplitude differs from the 3 Po~ m °del amplitude by a factor of 
l + fb/2/3 2 for both S + S and P + S final states. This explains the systematics of an earlier 
calculation || table II : compare columns 2 and 3] by Isgur and Kokoski, with j3 = 0.4 
throughout. Hence in the approximation where (the final state) /3 is constant throughout 
for all mesons, the flux-tube and 3 Po-models yield identical amplitudes when normalized to 
experiment. 

The reason why Isgur and Kokoski's results for the 3 Pq and flux-tube models are slightly 
different, is because they use a cigar-shaped overlap, which modifies the infinite cylinder 
overlap of eqn. |^ longitudinally away from the QQ-axis. The small deviation of their 
results for the two models from one another (~ 5% in width) demonstrates that this overlap, 
although physically reasonable, makes little quantitative difference on the scale of model 
errors, and is hence safely neglected, vindicating our choice of flux-tube overlap in eqn. |3]. 
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B Appendix : Decay amplitudes 



B.l Lowest radial states 

The lowest radial state wave functions were listed in eqn. |6|. Defining 

A- (^-A° M) 1 S;r 3/4 ^ exof ( M ) P * ) 

l 9^2 yioo V J (1 + /6/(2/? 2 ) ) (2f3l + pf 2 H [ ™ + M } 2fa + fP ) 

(11) 

we introduce the partial wave S, P, D and F definitions 

m + M M 



Si = V^^(6/3 2 (2^^ 



15 (2(3\ + (3 2 f m + M F " y r k r A rj ^ m + M FBiyH M 
F2 = -i-%=A ^ 5 2(— ^— pu) 3 (/? 2 + ^(2/?i + /3 2 )) (12) 

in terms of which we list the amplitudes in Table [5], consistent with the literature ||. 
For D** final states we employ a 3 Pi/ 1 P± mixing angle 9 according to the convention 

( l +L \_( cos0 sin^ \ / ^ \ ( . 

\ l +H ) ~ { -sine cos9 ){ 3 P 1 ) [l6) 

where 1 +L and 1 +H indicate the low and high mass states respectively. 
B.2 Radial excitations 

Instead of deriving the decay amplitudes for radially excited mesons afresh, it is more 
immediate to notice that the radially excited wave functions can be related to the lowest 
radial wave functions by differentiation |32j : 



^ {16 - 20/JV + 4,V } e-^ = -|{3 + 8^ + 8/^ 
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Table 5: Partial wave amplitudes M^(A — > BC) indicated in terms of the functions defined 
in eqn. [12] and named in accordance with partial waves L = S, P, D or F. Flavour can be 
incorporated for various final states by multiplying M L (A — > BC) by 1 (for D S D S , D*D*), 
y/2 (for DD,D*D*,D*D a ) or 2 (for D*D). 





M L (A -> BC) 






3 £l- 


* So 


1 S'o 




3 s 1 - 


> 3 P 2 


1 S'o 




3 ^- 


- 3 s, 


1 S'o 


Pi 


3 s 1 - 


» 3 Pi 


% 




3 S1- 


- 3 5x 


3 ^i 










Bx/V6 


3 P 2 - 


■> 1 S'o 


'So 




3 s l - 


» ^ 


% 




3 P 2 - 


- 3 s 1 


1 Sq 


D 2 /V2 










3 Bl- 


- 3 ^ 


% 


S2/V3 


3 A- 


-> 1 S'o 


1 S'o 


-B 2 








d 2 /Vq 


3 D 1 - 


- 3 5x 


1 S'o 




3 Po- 


■> 1 S'o 


x So 


-S2/V6 


3 D,- 


- 3 5 x 


1 S'o 


V2P2/V5 




- 3 ^ 


x So 


S2/V6 
















-D 2 /V3 











*p(r) = 4^/|^ r {5 - 20V} W« = 2^|/J"^,(r) 

* D (r) = ^ r 2 {7 - 2/?V} F 2Ml (f) e"^ = 2^A to(r) 



1 /? 3 / 2 1 9 9 

^ 45 (r) = - ^_^ 3/4 {105 - 210/3V + 84/3V - 8/3 6 r 6 } r 

" {2^ 2 3i + 24/5 4 -^ + 8/? 6 -^}^(r) (14) 



v/7! d 2 /? 2 ^ d 3 p 2 ' 

The differential operator depends only on /3, and can hence be pulled out when the 
integration in eqn. |2| over and y is performed. The amplitudes for radially excited 
meson decay can hence be found by applying the differential operators in eqn. |TJ] to the 
amplitudes in table ^|. 
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Figure Captions 

Figure 1: The contours indicate 0(3770) total hadronic widths fitting the mean experimen- 
tal width in eqn. |l|. Where there is 1D-2S mixing with mixing angle 8, we display sin# 
(which can be positive or negative). Pure ID and 2S states are also shown. The "square" 
indicates theoretically acceptable parameters. 

Figure 2: Total hadronic widths of the higher mass charmonia 0(4040), 0>(416O) and 
■0(4415) fitting experiment. For the latter two states the contours enclose narrow regions 
with experimentally acceptable total widths, as indicated by shading in one case. They in- 
dicate lcr variations from the mean experimental width (eqn. |l|). We display 0>(416O) as 3S 
or 2D, and -0(4415) as 4S. The regions in table |]where 0(4040) fits experiment are indicated 
by three large "dots" . The "square" indicates theoretically acceptable parameters. 

Figure 3: Total hadronic widths of the 2P states with J PC = 2++ (3.98 GeV), 1 ++ (3.95 
GeV) and ++ (3.92 GeV) |TB[ for theoretically acceptable parameters. We indicate regions 
where these states may be narrow. The lighter and darker shading indicate 2 ++ widths less 
than 5 MeV and 1 MeV respectively. The contours for ++ and 1 ++ enclose narrow regions 
where widths are less than 5 MeV or 1 MeV. 
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